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PREFACE 


This  paper  fxirnishes  an  explicit  mathematical 
treatment  of  the  Impact  forces  on  a  projectile  with  a 
spherical  nose  entering  water  vertically.   On  the  basis 
of  the  results  obtained,  a  less  technical  meraoranduin 
will  follow  stressing  the  more  Intuitive  physical  aspects 
of  the  problem.   In  addition,  an  extension  of  the  theory 
to  the  problem  of  oblique  entry  Is  being  considered* 

This  study,  related  to  v/ork  In  which  the  Applied 
Mathematics  Group  at  New  York  University  has  been  engaged 
at  the  request  of  the  David  Taylor  Model  Basin  and  of  the 
Naval  Bureau  of  Ordnance,  grew  out  of  Informal  suggestions 
made  by  Ensign  David  Gllbarg  of  the  Naval  Ordnance  Labora- 
tory.  The  work  was  carried  out  by  Dr.  Max  Shiffman  and 
Dr.  Donald  G.  Spencer  of  the  New  York  University  Group  of 
the  Applied  Ifethematics  Panel  with  the  assistance  of 
Mr*   H.  Goldstein.   Expert  numerical  help  was  rendered  by 
the  Mathematical  Tables  Project* 


R.  Courant 
Contractor's  Technical  Representative 
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I«  Introduction  and  Smmnary 

1#  The  Impact  force*   The  entry  of  a  bomb,  torpedo,  or 
mine  from  air  into  water  takes  place  in  three  stages: 
(i)  the  initial  stage  of  short  diiration  following  contact 
with  the  svirface  and  ending  with  formation  of  a  cavity  be- 
hind the  projectile;  (ii)  motion  through  the  water  with  a 
trailing  cavity  extending  to  the  surface;  (lii)  collapse 
of  the  cavity  and  sinking,  a  stage  in  which  the  cavity  is  a 
(possibly  pulsating)  air  bubble  following  the  projectile. 
Since  the  deceleration  experienced  by  the  bomb  is  generally 
largest  during  the  impact  stage  (i),  it  is  iroportant  to  ob- 
tain a  quantitative  estii;iate  of  the  impact  force  for  the 
purpose  of  determining  possible  damage  to  the  mechanism. 
In  this  report  we  are  concerned  solely  with  the  first  phase 
of  the  entry. 

Purthernore,  we  suppose  that  the  nose  of  the  projectile 
is  spherical  in  shape  and  that  the  initial  inrpact  with  the 
water  surface  is  vertical. 

The  upward  force  P  acting  on  the  sphere  in  the  impact 

stage  can  be  expressed  in  terms  of  a  dlmensionless  impact-drag 

coefficient  C  .   Let  A  be  the  radius  of  the  sphere,   U 
p  -^    '    o 

its  velocity  on  impact  (supposed  vertical),  and  p  the  density 
of  water.  Then  the  Impact  force  P  is  conveniently  expressed 
as 

P  =^pu2.TrA2.Cp 


where   C   is  the  "impact-drag  coefficient".   The  coefficient 

ir 

C     depends   on  the    time   -  It   is    zero  at   the   Instant  of   impact. 
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rapidly  increases  to  a  maximum  value,  and  then  decreases. 
It  is  the  main  object  of  this  report  to  determine  C  ,  and 
in  particular   the  maximvm  value   of   C  .      The   result   is: 

maximiom  value   of   C     =  •95. 

The  maximum  occurs  when  the   lowest  point   of   the   sphere  is 
about  one   quarter  of  a  radius  past   the   original  level  of 
the  water   svirface.      At  any  time      t      let     B     be   the   distance 
from   the   lov/est  point   of   the   sphere   to   the    original   level 
of   the  water    (see   Fig.    1).      The   graph  of      C        as   a  function 

IT 

of     B/A      is    shov/n  in  plate   2.      (Plates   are  placed  at   the 
end  of   the  report.  ) 


I-nitial  leoeC 


Figure   1 


Behind  and  alongside  a  rapidly  moving  projectile 
entering  the  water,  a  cavity  is  soon  formed.   It  is  knov/n 

for  the  case  of  a  sphere  that  this  cavity  is  formed  when 

B  T 

X  is  somev/here  between  -r   and  1.   Thus  the  graph  of  plate 

2  should  certainly  not  be  used  for  x  larger  than  1,  i.e. 

when  more  than  half  the  sphere  is  submerged.   It  may  be 
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pointed  out,  hov7ever,  that  plate  2  is  valid  for  the  whole 
range  of  -r-  if  no  cavity  is  formed,  which  occurs  if  the 
Initial  velocity  of  the  sphere  is  small. 

From  the  standpoint  of  damage  it  is  also  of  Interest 
to  know  the  average  pressure  p  acting  on  the  portion  of 
the  sphere  in  contact  with  the  water.   At  any  time   t   let 
a  be  the  radius  of  the  circle  of  intersection  of  the  sphere 
with  the  Initial  level  of  the  water  surface  (see  Pig.  1). 
Then 

P?  B 

where  K  =  •— r-  C  .   The  dependence  of  K   on  x  is  shown 
p   ^2  p  p     A 

in  plate  3«   It  is  noteworthy  that  the  value  of  K   at  the 

initial  contact  of  the  sphere  with  the  wab'or  is  infinite. 

(Actually  the  coiqDressibility  of  the  water  eases  the  entry 

somewhat,  and  the  initial  pressure  p  is  pcU   instead  of 

infinity,  where  c  is  the  speed  of  sound  in  water.) 

2.  Approxiraation  by  flow  about  a  lens.   Consider  the  flow 
produced  in  the  water  as  a  result  of  the  entry  of  the  sphere. 
As  the  sphere  penetrates  the  surface,  the  free  sui'face  no 
longer  remains  plane.   The  water  rises  slightly  and  spurts 
are  fortaed  running  up  the  sides  of  the  sphere.   The  exact 
shape  of  the  surface  is  difficult  to  deterinine;  estimates 
of  the  shape  and  its  influence  on  the  impact  coefficient  will 
be  presented  in  a  separate  memorandum.   In  the  present  report 
we  shall  neglect  the  motion  of  the  surface  and  likewise  neg- 
lect squares  of  velocities  on  the  surface.   The  boundary  condi- 
tion on  the  free  surface  is  then  0=0,  where  0   is  the 
velocity  potential  describing  the  flow.   The  potential  flow 
can  then  be  extended  to  the  infinite  space  on  both  sides  of 
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the   free   surface,   and  the   flow  at  each  Instant  corresponds 
to   the   flow  about  a   synnetrical   lens   as   Indicated  in  Pig.    2i 
It   is    the   flow  about  a   lens   which  fonns    the  main  matheraatl- 
cal  contribution  of   this  report.      It  v/ill  be  determined  by 
explicit  mathematical  formulae.      The  method  is  based  on  one 
used  by  A.    Som-nerfeld  in  the   theory  of   diffraction  of   light 
and  might  be  useful  in  the   analysis   of   other  hydrodynaralcal 
and  sound  phenomena. 


^  =  o 


Figure  2 


As  said  before,  the  flow  about  the  lens  provides  only 
a  first  approximation  to  the  actual  flow.   This  approximation 
is  poor  near  the  edge  of  the  lens  because  the  water  streams 
along  the  back  surface  of  the  lens.   In  the  actual  flow  the 
water  cannot  penetrate  the  remainder  of  the  sphere  and  rises 
instead  along  its  sides. 
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3«   Exact  foriaulat5.on  on  the  basis  of  the  momentum  equation. 

Let  B  =  D(t)  be  the  distance  at  time   t  of  the  lowest 
point  of  the  sphere  from  the  Initial  surface  level.   Let 

B0(x,3r,zjt) 


Figure  3 


be  the  velocity  potential  describing  the  flow  at  the  Instant 
t«   Here  B  =  -^   ,  and  x,y,z  are  rectangular  coordinates, 
the  2-axls  being  vertical  and  z  <  0  being  the  half  space 
initially  occupied  by  the  water.   The  vertically  downward 
moraentUM  of  the  water  is 


(3.1) 


m 


pB|||r^  dxdydz   =     pB     0dxdy     -  pBL, 

W+S 


^S* 


Here  W  Is  the  wetted  portion  of  the  sphere,  S  Is  the  free 
surface,  and  L  Is  defined  to  be  the  limit  of  the  integral 
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//• 


0dxdy 


extended  over  the  submerged  portion  H  of  a  large  spherical 
membrane  the  radius  of  which  tends  to  infinity. 

The  external  forces  acting  on  the  fluid  contained  in 
the  membrane  H  are:  (i)  the  downward  force  of  magnitude 
P  exerted  by  the  sphere;  and  (ii)  the  upward  force  exerted 
on  H  by  the  pressvire   of  the  fluid  external  to  it.   (The 
resultant  force  due  to  the  constant  atmospheric  pressure  is 
zero.   The  only  effect  of  gravity  is  the  buoyancy, which  is 
negligible  compared  to  the  impact  force  P.)   The  force  (ii) 
is  equal  to  the  time  rate  of  change  of  the  term  -  pBL  in 
(3.1),  as  may  be  seen  from  Bernoulli's  equation.   Since  the 
resultant  of  the  external  forces  is  equal  to  the  time  rate 
of  change  of  the  momentum  (3. 1),  we  therefore  have 

(3.2)  ^  =  -^  (^^) 

where 


(3.3)  M  =  p 


0dxdy. 


W+S 


The  quantity  M  may  be  considered  to  be  the  effective  "momen- 
tum mass"  of  the  water. 

Let  M^  be  the  mass  of  the  incoming  projectile.   Since 

(3.4)  -P=^(M^B)  , 
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we  obtain  from  (3»2) 

(3.5)  (M  +  M^)B  =  M  U^    ,    or  B  =  -^ 

o  o  o    '  2.+fi 


where      /i  =    •«•  •      Equation   (3.5)    expresses   conservation  of 
moraentura. 

Prom  (3.4)   and  (3.5),    it   follows    that 


dM  di'I 

(3.6)  p   =  -dB.  J      =       ^    ,U^   . 

1+//  (l+/i)^     ° 


Introducing  the  dirnenslonless   quantities 


(3.7)  m  =-JL_   ,      b   =  I   , 

|vrpA^ 


we  may  write 


(3.8)  P   =  ^PU^-TTA^.C 


xvhere 

dm 

(3.9)  C     =        ^ 


For  Incoming  projectiles   of   large  mass   compared  to   the 
mass   of  water  displaced  by  the  hemispherical  nose,    fi    may  be 
neglected.      Then  the   impact-drag  coefficient     C       is   merely  ^ 


CONFIDENTIAL 


-  8  -  CONFIDENTIAL 


The  problem  Is  thus  to  find  M  or,  since  the  exact 
flow  l3  difficult  to  determine,  a  reasonable  approximation 
to  M. 

The  result  obtained  in  the  present  paper  is:  the 
dimensionless  virtual  mass  m  and  the  impact-drag  coefficient 
^     depend  on  the  depth  of  penetration  b  in  the  way  indica- 
ted in  plates  1  and  2.   The  maximum  value  of  -gr  is  .95,  and 
it  occurs  for  b  =  •24. 

It  is  possible  to  use  an  energy  equation  in  place  of 
a  momentum  equation,  but  this  would  require  a  more  detailed 
analysis  of  the  free  surface.   For  the  spurts  that  run  up 
the  side  of  the  sphere  may  carry  considerable  energy,  whereas 
the  vertical  momentum  they  carry  is  small.   Also,  there  may 
be  energy  losses  due  to  sound  waves  in  the  water  and  elastic 
vibrations  of  the  projectile. 


4.   Conqparison  of  lens  flow  with  disc  and  ellipsoid  flows. 
Various  writers  (see  references  [2],  [4]  and  [5]  in  the  bibli- 
ography at  the  end  of  the  report)  calculate  C   on  the  basis 
of  the  following  assuinption:  the  entering  body  at  each  Instant 
of  time  t  is  replaced  by  its  cross  section  at  the  surface 
(assumed  to  be  the  initial  surface)  and  the  virtual  mass  M 
is  taken  to  be  that  due  to  the  flow  about  this  disc   If  the 
disc  has  a  radius  a,  the  virtual  mass  Is,  from  classical  hydro- 
dynamics, equal  to 

4   3 
•^pa  . 


The  dimensionless  virtual  mass  m^^  is  then 


«^=i^=^[b(2.to))3/2. 


is'" 
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and  the  impact-drag  coefficient 


dm. 


db" 


is 


^   =|(l-b)V'b(2-b) 
do    Tf 


A  comparison  of  this  approximation  with  the  lens  flow 

approximation  used  in  this  paper  is  given  in  plate  4,   The 

dmi 
agreement  is  fair.   The  maximum  of  -^^   in  the  disc  approxi- 
mation is  1.27  compared  to  .95  for  the  more  accurate  lens. 

Instead  of  a  disc,  consider  an  ellipsoid  of  revolu- 
tion which  fits  the  lens  as  shown  in  the  diagram,  with  semi- 
axes  a,  a,  B.   The  flow  about  such  an  ellipsoid  would  seem 
to  be  a  better  approximation. 

The  flow  about  an  ellipsoid  is  classical  (see  [l], 

pages  456  and  457).   Write  y  =  —  and  assume  that  0  <  y  <  1. 

a 


Z,<£'>7v3- 


Sll  tpsoiJ  ^ 


Pigtire  4 
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The   dlnensionless  virtual  mass     nig      then  is 


2   -  a 
o 


where 


«o  = 


=:  — =— -  (  1 arc  cosy  )       (0  <  arc  cos  y5  i ) 

and  the  impact-drag   coefficient   is 

^  -  .  4^.,.    _    ..       «o     ^  8  1/^7^-^   {       y      \^  Ifo 
-^-  -3b(3b-  4)   ^— -  +  3  Vb(2-b)   ^^-^j     —     . 

It  is  seen  from  plate  5  that  the  agreement  with  the  lens  flow 

diTip 
Is  better   than  that  of   the   disc.      The  maxinum  of  -rrr-  in  the 

ellipsoid  approximation  is    1»1,    compared  to   ,95  for   the    lens. 

As   another  possible   approximation,    take   a   sphere  with 

center  at   the   center  of  the   lens   and  radius   equal   to     a.      The 

1  3 

virtual  mass   is     -t^PR    »    and  the   impact-drag  coefficient  is 

J  =  (2-2b)(2b-b^)V2  . 


■crr 


dmj 


A  graph  of  -w-  is  given  in  plate  6,  and  it  is  seen  to  be  a 
better  approximation  in  the  neighborhood  of  the  maximum  than 
the  corresponding  quantity  for  either  the  disc  or  the  ellipsoid. 

We  have  been  concerned  in  this  paper  only  with  vertical 
entry.   However,  for  oblique  entry  it  would  seem  desirable  to 
replace  the  lens  by  an  ellipsoid  and  then  to  solve  this  somewhat 
easier  problem. 
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II»    Mathematical  Treatment* 

5#      Preliminary  remarks*      The  f lov/  pattern  of   the   fluid  at 
the   time      t     is    described  by  a  velocity  potential 

B0(x,y,zjt) 
which  satisfies   the  Laplace   equation 


d>  „      +<6         +0         =0 

^xx       ^yy       ^zz 


Figure  5 
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together  with  the  following  boundary  conditions: 

(a)  On  the    sphere,    the   negative   of   the   normal  deriva- 
tlve   of     B0  is   equal  to   the   cojiiponent  of    the  velocity  of   the 
sphere  along  the  normal,    or 

-  B  1^  =  B  cosX. 
8n 

(b)  By  virtue  of  the  assuriiptions  made  in  Part  I,  Sect.  2, 
0  vanishes  on  the  plane   z  =  0« 

To  see  how  to  construct  this  potential  function  0  ,  con- 
sider first  the  case  in  which  the  sphere  is  completely  imnersed 
in  the  'vater  and  moving  vertically  (see  Fig.  5),  A  classical 
method  of  constructing  the  desired  pot^-^ntial  function  is  "by 


-11  =  ^''^^ 

OTi 


Figure   6 
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the  method  of  Images  (see  [1],  page  420).   Suppose  that  the 
free  siirface  were  not  present.   The  flow  produced  by  a  sphere 
moving  in  an  infinite  body  of  fluid  is  that  created  by  a 
dipole  of  suitable  strength  placed  at  the  center  0  of  the 
sphere.   This  satisfies  the  boundary  condition  (a)  on  the 
sphere.   However,  because  of  the  presence  of  the  free  surface, 
the  boundary  condition  (b)  on  the  plane   z  =  0  is  violated. 
To  satisfy  (b)  introduce  the  image  of  the  dipole  relative  to 
the  plane   z  =  0,  which  is  a  dipole  of  equal  strength  at  the 
reflected  point  Q,.   The  condition  (a)  on  the  sphere  is  now 
violated,  and  to  remedy  it  Introduce  the  image  of  the  latter 
dipole  relative  to  the  sphere.   This  image  is  a  dipole  of 
suitable  strength  at  the  inverse  point  Qg*   Indefinite  con- 
tinuation of  this  process  leads  to  a  potential  function  which 
satisfies  both  the  boundary  conditions  (a)  and  (b). 

We  shall  apply  a  generalization  of  the  above  process 
to  the  case  where  the  sphere  intersects  the  free  s\irface  (see 
Pig.  5).   The  process  as  it  stands  leads  to  a  serious  difficulty, 
namely  that  some  of  the  successive  Images  lie  in  the  fluid  be- 
low the  sphere,  thus  creating  dlpoles  (singularities)  inside 
the  fluid  which  violate  the  assumptions  of  the  problem.   For 
example,  in  Pig.  5  (where  B  <  -s),  the  image  of  the  dipole  at 
the  center  of  the  sphere  relative  to  the  plane  z  =  0  introduces 
a  dipole  in  the  fluid. 

The  difficulty  may  be  overcome  by  using  multl -valued 
dlpoles,  the  Images  of  which  lie  in  the  various  replicas  ("sheets") 
of  a  Rleraann  space  and  do  not  give  rise  to  singularities  in 
that  portion  of  the  space  in  which  we  are  Interested.   The  con- 
struction of  the  multi-valued  dlpoles  is  based  on  the  method 
used  by  A.  Sommerfeld  in  his  classical  work  on  the  theory  of 
diffraction  (see  [3]). 

We  describe  the  Riemann  space  in  the  next  section,  and 
discuss  the  multi-valued  dlpoles  in  sections  7-9.   The  potential 
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function  $  satisfying  the  boundary  conditions  (a)  and 
(b)  is  constructed  in  section  10,  the  final  formula  being 
(10.11). 

The  velocity  potential  0  ,  extended  to  the  whole 
space  by  reflection  in  the  plane   z  =  0,  describes  the  flow 
about  a  syinmetrical  lens.   In  Appendix  A  we  discuss  the 
flow  about  an  asymmetrical  lens  formed  by  the  intersection 
of  any  two  spheres. 

Finally,  if  the  sphere  is  entirely  below  the  surface 
z  =  0  and  does  not  intersect  it,  the  construction  of  0  is 
classical  (as  stated  above).   In  Appendix  B,  however,  it  is 
noted  that  the  potential  function  can  also  be  obtained  merely 
by  replacing  some  of  the  real  quantities  occurring  in  our 
formulas  by  purely  Imaginary  ones. 

6»   The  coordinate  system  and  the  Riemann  space.   We  select 
a  coordinate  system  which  fits  the  problem  better  than  either 
the  ordinary  rectangular  or  cylindrical  system.   It  is  ob- 
tained from  bipolar  coordinates  in  the  plane  (see  [l],  page 
165)  by  rotation  about  an  axis. 

Let  r,z  be  ordinary  rectangular  coordinates  in  the 
plane.   Introduce  in  place  of  r,z   the  bipolar  coordinates 
a  ,  ip   defined,  in  complex  terminology,  as  follov;s: 

(6.1)  a  -  1  ^  =  log  Il-jL-LUl-S.  ,   or  r  +  iz  =  al  cot  i(  lA  +  lo). 

r  +  iz  -  a  '^ 

Geometrically,  a    is  the  logarithm  of  the  ratio  of  the  distances 
from  the  two  points   -a  and  a  on  the  real  axis  and  ip    is 
the  angle  subtended  by  the  line  Joining  these  two  points  (see 
Fig.  7). 

We  observe  that  the  curve  ip  =   constant  is  a  circular 
arc  connecting  the  points  P,  and  tg.   The  curve  a    =  constant 
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is  a  circle  with  center  on  the  real  axis  and  with  the  points 
P,  and  Pg  as  mutually  inverse  points. 

The  value  of  ijj   on  the  upper  side  of  the  line  segment 
P^Pp  is  tr  ;  the  value  of  i^  on  the  lower  side  is   -  tt. 
The  simple  plane  is  obtained  by  identifying  these  two  values 
of  the  angle  i/*  .  But  if  the  angle  \p   is  allowed  to  vary 
outside  the  range  -  v    <  \p  <  v  ,     and  if  values  of  0  differ- 
ing by  2iT  are  not  identified,  we  obtain  a  Riemann  surface 
with  P,   and  Pg  as  branch  points.   The  Riemann  surface 
which  we  consider  is  the  infinite-sheeted  one  in  which  i//  is 
permitted  to  take  any  value  from  -oo   to  +00  with  no  identi- 
fications. 


z 

(r=  0 

^  P 

-a -- 

n^=  -h'Tf          a.  / 

r 

^ 

(T-    Uj>    ^^' 

Pigxu^e   7 


CONFIDENTIAL 


-  16  -  COIIFIDENTIAL 


A  system  of  coordinates  in  space  is  obtained  by 
rotating  the  plane  about  the  z-axis«   Let  0    be  the  angle 
of  rotation.   The  space  coordinates  (^  ,  }p  ,  6   are  called 
toroidal  (because  the  surface  a   =  constant  is  a  torus). 

To  the  Riemann  surface  discussed  above  there  corre- 
sponds in  space  an  infinite-fold  Riemann  space.   The  branch 
points  at  ±  a   are  replaced  in  space  by  the  circular  branch 
curve  r  =  a,   z  =  0,  and  the  cut  line  joining  P,  and  Pg 
is  replaced  by  a  cut  svu'face  r  <  a,   z  =  0  across  which  the 
replicas  (or  "sheets")  of  the  space  are  joined. 

The  equation  (6.1)  defining  bipolar  coordinates  can 
be  expressed  in  real  fona  as  follov/s: 


(6.2)       r  =  a  ^^^^^^^  ,      z  =  a  ^l^Jl 

3  -  r    '  3  -  T 


where  s  =  cosh  a  ,   t=  cos  i^   .   The  element  of  length 
in  space  is 

2 

(6.3)  dz^  +  dr^  +  r^d6)'^  =  — ?,  (da^+  di/;^  +  slnh^a  dO^)  . 

(3-t)2 


We   observe  finally  tlmt   the  points    {  a  ,  ip    ,  0  )   and 
(  ^  >'/'o»  ^  ^    ^^^   inverse  with  respect   to   the   spherical   cap 
(p    =  constant  =  c     provided  that 

(6.4)  rp^  +  xjj^  =  2c. 

For  the  mapping  in  which  the  point  {<^  ,ip)    goes  into  (  <^,^+  c), 
where  c  is  a  fixed  real  quantity,  is  a  linear  transformation 
of  any  plane  9   =  constant  Into  itself.   Since  the  above  state- 
ment about  inverse  points  is  true  when  c  =  0,  it  is  true 
generally  because  it  is  Invariant  to  the  linear  transformations. 
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7.      Dlpoles*      In  this    section  we   obtain  the   expression 
for  a   dlpole    in   toroidal   coordinates. 

The  Inverse   distance  between  a  point  on  the   z-axls 
<T     =0     at     i//  =  lA-      and  an  arbitrary  point    (  <t,iI^  ,9)    of  the 


space   is,    froin   (6.1) 


1  1  |sinW°l  /  g       _     -. 

(7.1)^  ^  ^  ./  ST 


|(r  +  iz)-(iZQ)|  a  Ks-   cos(  il/  -  ip^) 


This  represents   a  unit  Blpole  on  the    z-axis    (cr=  0)   at   the 
point     z  =  z     =  a   cot  ^'-*  •       Instead  of    (7.1)   we   introduce 
the  function 


(^.-S)  P^      =  3in|o|/ 3_Z.jr 

^O  2     ^   S     -     C0S(I//  -    i/;^) 


which  represents   a  pole   of  strength  la. 

A  dipole  at   the   same  point   is   obtained  by  differen- 
tiating   (7»2)   with  respect   to      z   .      However  we   shall  find 
it  more   convenient   to  use  a   dipole   of  suitable   strength, 
namely 

(7.3)   qAa,ip)  =   -  ^     esc     -     ^   (p^^)    =  a   esc  ^  ^   (p^^) 
^o  o  ^o 


a 


1^ 
(s-T)     sln(^-.A)       ^^^  ^o  (3    -r)^^^ 

[s-cos(r/;-  lAo)]^  ^     [s-cos(^-  ^J]^^^ 


This  function  is  a  dipole  of  strength  -  §  |csc"^  ^°I  situated 
at  the  point   z  =  z^  =  a  cot  ^  on  the  z-axis. 
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The  Stokes'    stream  fu.nction     tjh        corpesponding   to 
^ipQ      is    (see   [  1]  ,    Chapter  XV,    especially  page   414): 

iry    A\      4-,      -  a^lo<=^3    ^oi    r^   _   1  ^2  siiih^  a 

(7,4)       t^     =  -    lose     w        -rr  =  w  a      ^Tp 77^    • 


The  functions  q^       and  tw,   are  connected  by  the  linear 
differential  equations 

fVo  ^  _  1  .^ 
da  r   3^ 

(7.5) 


3)//      r   3  <T 


Y/e  wish   to  find  the   image   of   the   dlpole      q^^     in  the 
spherical   cap   ip  =   constant  =   c,      that    is,    vie  wish  to  find 
a  potential  function    0    such  that    the   normal   derivative    of 
q^     +   0  is   zero   on    0  =  c.      This    condition  is    equivalent   to 
the   stafeenent   that    the    sum  of   the   corresponding  Stokes' 
stream  functions   is    constant   over   the   cap.      We   sliall   show 
that   the   image    is      -  q^'      where     i/^     +  j/;'    =  2c,    or   c   -     (/'q  = 
-   (c   -!/»')•      The  functions      q^/.       and     q^J      are   dipoles 
situated  on  the   z-axis   at   the  points   given  by    ^  -  4^^  and 

ip  =  ipl^    ,    and  these  points   are  mutually  inverse  relative    to 
the   spherical   cap    ^  =   c.    (See    (6«4)») 

The   Stokes'    stream  function  for    the    sum     q,/,     -   q^ 
is      t^     -   t^,'      which  vanishes   on    0  =  c  by   (7.4)    since 
cos(c   ~  ^q)   =  cos(c    -  ^i)»      Thus    the   normal   derivative   of 
qw,     -  qj/,'      is    zero   on   i//  =   c»      (This    is    the   reason  for   taking 
a   dipole   of   strength     -  ^   |  esc'  ^°  |    as    the  basic   dipole      qw,^. 
The   stream  function  then  has   a   simple   form  and   leads    to  a 
simple   image   result.) 
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8«  Multi-valued  dlpoles»  The  construction  of  our  multi- 
valued dipolea  will  be  based  on  the  method  of  Sorainerfeld 
referred  to  above.  Since  ip  is  involved  in  an  analytic 
way  in  Pj/.  and  (l^ij  ,  these  functions  of  a  and  i//  still 
satisfy  the  potential  equation  when  ip  is  replaced  by  a 
complex  nxamber  ^  •  They  may  be  regarded  as  representing  a 
pole   or   dipole   respectively,    situated  at   the   "complex  point" 

a  =  0,    iA  =  I    •      Also      tt      Is    the    stream  function  correspond- 
ing to      q.      ,    since   the   differential  equations    (7.5)    are   still 
satisfied  when     ^    is    complex* 

As   functions   of     ^    in  the   complex     ^ -plane,      p^    ,q.    ,t* 
are  multi-valued  and  have  branch  points   at   the   zeros   of 
s-cos(  i/"  ^  ) ,  that    is   at    |  =  ip±  la  +  2kTr  where  k=0,  ±1,  ±  2, 


.  • .    • 


i    Y-S'Tf  ■i-i<r 


S  - piive. 


6 


w 


"y^  +  L<r 


If  ■y/'-t-^TT  -t-ixr 


r2\    '^-ir 


W 


Figure   8 


CONFIDENTIAL 


-   20   -  CONFIDENTIAL 

Cut    the     I -plane  by  vertical   lines  proceeding  from  the  hranch 
points   to    oo    as   shown  in  the   diagram  (Pig.  8).      The   functions 
Vf.       f<li      f    tt        are   now  single-valued  in  this    cut       | -plane 

and  take   opposite   signs   on  the   opposite    sides   of   each  cut.   V/e 

1/2 

define    [3-cos(!/'-^  )]    '       to  be  positive  v/hen    ^    is   real. 

The   basic  multi-valued   dipole   situated  at  a    (real)   point 
,    ip  =  ip        on  the   z-axi; 
space    is    defined  as  follows: 


c    =  0,    0  =  ^        on  the   z-axis    of   our   infinite-sheeted  Riemann 


W 

where   the  path     V/     of   integration  is   the   two-branched  path 

indicated  in  the  diagram.      The   function     q,,      (<y,^)    is   a 

potential  function  since   differentiations    can  be  performed 

under   the    integral  sign,    the  path     W     being  teroporarily  fixed. 

(For  an  n-sheeted  Riemann  space  in  which  values    of   the 

angle    i//    differing  by     2nTr       would  be   identified,    the   corre- 

sponding  n-valued  dipole   q^.   (o-,?^)  is 

^o 

(8.2)  q(n)^(a»  =^Jq^  .^cot4^°d|   . 

w 

^  -^o 
Since  cot  — is  periodic  in  ^  of  period  2nTr  ,  and  q^ 

is  periodic  of  period  2v    ,  it  is  easy  to  see  that   q,  (^,^^) 

is  periodic  in  ip    of  period  2n  tr  .  )  ° 

We  now  briefly  describe  the  properties  of   qv^('',^)• 

(A   similar   description  applies    to      q!l(<'',^)»)      On  acc°ount   of 

the  pole   of   the    integrand  in   (8.1)    at  °  ip       one   sees,   by  the 

Cauchy  residue    theorem,    that 

(8.3)  (i^^-     q^^(a,0) 
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(oo) 

Is  regular  as  a  ^  0  and  ^ ->  'PA^q¥=2^v),      Thus   q  , 

behaves  like  a  dlpole  near  a  =  0,  ip  =   0  •  But  (f"?'  {Ofip) 
is  regular  near  a  =  0,  i// =  ^^  +  2kn,  k  =  i  1,  ±2,  .°  . ,  since 
the  Integrand  Is  regular  at  these  points. 

In  addition  there  are  still  the  poles  of  q^  Itself 
at   ^  =  2kir,  k  =  0,  ±  1,  ±2,  •••  due  to  the  presence  of  the 
terra  cot  4-  (see  (7»3))«   One  easily  sees  that 


as    (T   ^  0     and        ip  ->  2krr,    k  =  0,   ±   1,   ±  2,    .. .    (  ip^^2iv)    . 

Thus,    the   function     q^°|^'  is   not  simply  a  dlpole  at  ip       but 

^  o 
it  has   certain  non-vanishing  values   at  oo  (  a  =  0,  i/'=  2ktT ) 

in  the  various    sheets  and  has   no   other  singularities.      (A 

constant  value  at    oo    may  be    thought  of  as  a  pole  at    <»    •      If 

in  the   definition  of  q.         the   tern  cot  X  were  replaced  by 

lb  ? 

cot  »°,    the  resulting   function  in   (8.1)   would  be   zero  at   the 

points  at    oo    .    this  would  be   a  simple   dlpole). 

In  the  above  wo  have   inrplicitly  assumed  tliat     ip       is 

not  a  multiple   of     2tt,   for  then  the    integrand  in  (3.1)   has   a 

multiple  pole.      Consider,    for  example,      (l^^{(^,^)»      The   Integrand 

in  (8.1)   has  a  double  pole  at    ^  =  0,    and  the   residue   there   is 


a   sjbi^.  „ 
2   s-Tr-^^ 

Therefore 


a  /       3   -  T         '\l/2 
^(s-co3(  ip  -^  )j 


-  o    3  3.n^    _  , 

=  a -r  =  z 

s    -  T 


(8.5)  qo(<^,«/')    -   z  -^  0 


CooX 


as    a  ^  0       and      i/' ->  0,    and  so      (\AOf^)   behaves    like      z     near 
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the  point  at  oo    in  the  first  sheet.   A  similar  analysis 

applies  to   ^^J^^^'^^* 

In  the  n-3heeted  case,  it  is  of  interest  to  consider 
the  special  case  n  =  1  where  we  are  in  the  original  single- 
sheeted  space.   The  integrand  in  (8«2)  is  of  period  2n   , 
and  a  simple  analysis  yields 

(8.6)  q^J(<r,0)  =  q  ,  {a,rP)    -  ^  cot  ^°  ^^^t  0  , 

so  q  ,  differs  from  an  ordinary  dipole  q  ,  by  a  constant. 
If  in  (8»6)  we  now  let   i/;  ^  0,  we  obtain 

It  may  also  be  verified  directly  that 

(8.7)  qo^(<^,^)  =  z  . 

9.      Inversion  of  multi-valued  dipoles.      It   is  well-known  that 
in  ordinary   space   the   image   of  a   dipole  with  respect    to  a 
sphere    is    another   dipole   situated  at   the   inverse  point.      In 
this   section  the   theore]n  will  he   extended  to  multi -sheeted 
spaces • 

Theorem  9.1.      The   function     q*T'(<^,^^)    -  q'T/Co',^)    has 

vanishing  nortnal  derivative   on   the   spherical   cap    1//=  constant  =c 

Proof;  Because  the  relations  connecting  stream  function 
and  potential  function  are  linear,  the  streaia  function  t^T'C^r.^) 
corresponding  to      q  ,  (a,^)    is 


t,    (<r,0)   =  -5^     tfc      .   __L_-.d^ 
W 
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v/here      tt      Is    the    stream  function  corresponding   to      q.    , 
obtained  from   (7,4)   by  replacing    i//^     by    ^  .      It   is  re- 
quired  to   show  merely  that     ^jj   {a,ip)    -  t^'(o-,(//)    is    constant 
on     tp  =  c» 

It   is   more   convenient   to   introduce  a   new   complex 
variable    ^     defined  by       ^  -  i/;    =    ^  ♦      Then 


(9.1) 


a 


-   (3-T)^/^slri^    ^  ^-^iP        (s    -t)- 


V2 


(s-cos^") 


372 


2        (s-cos^)V2 


(9.2) 


t     -  1  a^ 
t^    -  ^  a 


slnh^tr 


(s-r)V2(s.cosa 


3/2 


and 
(9.3) 


qX(-,^) 


1 

2TTi 


W 


^   +    i//  -  ^^, 


d^ 


(9.4) 


w' 


d^ 


wheiB    W    is   the  path  in  the     ij" -plane   shown  in  the   figtire. 

I 


I 


I 


t 
I 


^ 


i<r 


r^-i<r 


'    ^' 


f    .  .'    .  ,' 


//=/^^;<£ 


Figure  9 
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By  the  symmetry  of  the  path  W',  the  lower  path  Wj_ 
is  obtained  from  the  upper  path  w'  by  replacing  ^  by  -f  • 
The  function  tt   is  by  (9»2)  an  even  function  of  ^,    There- 
fore (9»4)  can  be  written  in  the  form 

(9.5)   t7^(-,'A)  =  ^J  t^   '(  t^i^-rp^)    ■"    ^-(/-^,)h- 

Likewise  for   t,'  {<^,4f)»      Evaluating  at  i//  =  c   and  noting  that 

the  hypothesis  \b^   +  \b  =   2c  can  be  written  c-iA  =-(0-0), 

o     o  j^,      o   ^^,  ^    _:o'» 

equation  (9»5)  shows  immediately  that   t^  (<^,C)  =  t^/  (<t,  c)  , 
and  the  theorem  is  proved.  °  ° 

Theorem  9«2«  The  function  q   (a,i/»)  +  q  .  (  <t,^)  vanishes 

^  ,co)        '^o  "^o 

on  \p  =  0»      In  particular  (lAo fXp)   vanishes  on  ^  =  0. 

Proof;   For  ip  =  0   the  function   q^  is  an  odd  function 

of  ^  by  (9.1).   The  same  argument  as  in  the  proof  of  theorem 

9.1,  applied  to  (9.3),  yields  q*^'  (a,o)  =  -  q"\  {a.O),      The 

^o  ~^o 

second  statement  of  the  theorem  is  an  immediate  consequence  of 

the  first  one. 

10.   Construction  of  the  potential  function.   We  are  now  pre- 
pared to  apply  the  method  of  images  discussed  in  Section  5. 
The  desired  potential  (p    is  regular  at  °o  {a=  0,^  =   0)  and 
satisfies  the  following  boundary  conditions: 

(10.1)  0  =  0   on    ^  =  0 

(10.2)  -^=cosX=:-|2    on   «^  =  -  ;5  . 

on  dn  ^ 

The  potential  function  z   satisfies  (10.1)  and  (10.2) 
but  is  not  regular  at  c»  .   Set 

(10.3)  V  =  ^  -  2 
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^  ~  o 


^=-/3 


=  cas  ?>,  = 


drt 


Flgtire   10 


Then     V     is   a  potential  function  satisfying  the  following 
conditions : 


(10.4) 


V  =  0        on  rp  =  0 


(10,5) 
and 


-    rr-       =    0  on 

3n 


ip  =  ~  /3 


(10.6) 


V     beliaves    like      -  z     at    <»• 


(The  potential  functioa     V      corresponds    to  the   flow  about  a 
fixed  lens,    the  fluid  having  unit  velocity  at  =».) 
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We  shall  construct  V  by  the  method  of  images,  using 
the  infinite-sheeted  Rlemann  space  described  by  0  ^  a  <  oo  , 
-oo<'/'<+oo.   (It  may  be  pointed  out  that  if  the  angle 
^   -  "^     where  n,  q  are  integers,  it  would  be  sufficient 
to  use  an  n-sheeted  Riemann  space  in  place  of  the  infinite- 
sheeted  one.  ) 

(OO) 

To  construct  V,  first  start  with  -q  {a,i//)   which  has 
the  correct  behavior  (1O.6)  as  a->   0,  ^  ->  0  (by  (8.5),  and 
satisfies  (10.4).   To  satisfy  (10.5),  one  must  add  q^/a(ff,j/;) 
by  theorem  9.1j  to  satisfy  (10.4),  one  must  now  add 

Cool 

qgitr,^)  by  theorem  9.2;  etc.   The  result  is 

(10.7)  V(a,^)  =^  (-1)^*H;L(^,'/') 

where  the  sum  means   lim  ^^^1^        .   We  demonstrate  the  conversieaice 
of  this  sum.  below  by  evaluating  it.   We  may  write  (10.7)  in  the 
forms 

(10.8)  V(a,,/;)  =  -q~'(^,^)  +^(-1)^"'^  fc;S^^'^)  "^  ^-Sk^^^'^lJ 

and 

00        I—  _ 

(10.9)  Vias^)    =  2.^-l>     ^2^^'^^    -   q-2(k+l);S^^>H    * 

Equation  (10.8)  shows  that  (10.4)  and  (10.6)  are  satisfied; 
equation  (10.9)  that  (10.5)  is  satisfied. 
Substituting  f«m  (8.2)  we  obtain 


vK.)-^J,,  {^c-.-^T^ 


w 
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or 


(10.10)  V  =-^Jq^  .  I^csc^d^  =.^jq^.  2!^esc|^(r+Od^ 


W  V 

Expressed  as  a  real  integral 


1  .  sinfl  t*"  I       3  -  T   tPj.^     tt^  V         irk   a.Trtl 
^  =  ^^-V'jjcosh  U  >  3  ^[°^^  I  ^3°  ^  +  ^°«°  2?  °°*  ^  ] 

where  ^  =  i//"  +  lu«   Th©  desired  potential  function  is  $  =  V  +  z 
by  (10.3).   It  is  easily  verified  from  (8.7)  that 


and  so  w©  may  express  the  desired  potential  function  as  a  real 
integral  in  the  form 


du 


( 10.  11)  $  =  M^ry^^i^  .  in.  [|  cot  I  C3C  ^p*  ^=^  flcot^L 

du 


-  2   CSC     -w 


We  observe  that  the  function  whose  iinaginary  part  appears  in 
the  integrand  of  (10.11)  is  bounded  at  ^  =  0,  and  we  therefore 
see  directly  that  ^=0  atoo(<r=0,  -^^  =  0)  because  the 
factor  ys  -  T  vanishes. 

11.  Virtual  mass  of  the  fluid.  As  stated  in  Section  2  we 
neglect  the  motion  of  the  surface.   Then  $  =  0  on  the  surface 
S,  and  the  wetted  portion  W  of  the  sphere  is  the  spherical 
cap  '\|/'=  -^  .   By  (3.3)  the  virtual  mass  of  the  fluid  becomes 
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where      dS      denotes   an  element   of   area    on'^=   -(b       and   A  is 
the   angle   indicated  in  Fig.  10.      Since      cos    -^   =   -  g^  ""Sn      * 
we   observe  that  under   these   circumstances      M     is    also   the 
kinetic-energy  virtual  mass    of   the  fluid. 
Since  $  =  V   +   z     v/e  have 


(11.1) 


M  =  M^   -  M^ 


where 
(11.2) 

and 


p  0 


M^  =  e 


V   cosA    dS 


^.-^ 


0  n 


(11.3)      M^   = 


-  c 


z   dxdy  =  -TT^A      cos     1   (2   -   cos  (3) 


Ajf^-^ 


-iTr^A'^(2   +  3  cos^-   cos^f) 


by  an  elementary  calculation.      In  fact      M_,      is   equal   to   the 
mass    of   the   fluid  displaced  by  the   sphere. 

The    element   of  area      dS      on  the   spherical   cap  ->!/■=   -6 


is 


^o        o«„      a       ,_  _  o«o2  sinho-  do-     _  ^^2        d  s 

dS    =    2TTr d<r  =    2Tra      — -rr-   =    2TTa      -. -rry 

s  -  T  (g   _t)2  (s   -V'^ 


z  -  z. 


where  t-  =  cos  ^     .      Also     cos  X   =  -  ■^  = T —     where      z 

is    the   z-coordlnate   of   the   center   of   the    sphere,    and     z      that 


CONFIDENTIAL 


-  29   -  CONFIDENTIAL 

slnC 


of  a  point  on  the   s\irf  ace  of   the  sphere.      Thus      z  =  -  a  ^    _  ^ 
and     2     =  -  cot  fi      and  we  have 

C03  A  =    -   ||=-|  (cot  p  -  |%£-)  =  -    Smp  .  (cot  p  -  |%^)(0  <    p   ilT  )  . 

Hence 

(11.4)  M^  =  -  2Trea2j    v((r,  -  (i  )^^^-^(cot  ^  -  |i^   )d  s   . 

To  evalroate  this  integral  it  is  convenient  to  put  the 
integral  (10.10)  for  V  in  a  slightly  different  form.   For 
■y\r  =  -  ^      the  Integrand  in  the  second  integral  of  (10.10)  does 
not  have  ^  =  0  as  a  singularity.   The  path  W  can  therefore 
be  deformed  into  the  path  W"   shown  in  Fig.  11,  and 

(11.5)  v((r,-^)  =  ^[  q, .  jrS  ^«<^  ?|  d^ 

where  '\[^   is  set  equal  to  -p  in  q,   •   Substituting  from 

(11.5)  into  (11.4)  we  obtain  (the  inversion  of  the  order  of 
Integration  is  easily  justified) 

(11.6)  M,=  .2.ea^.  -1-  [(  r.^li^^(eotp.f^±)ds)^^^^ 

W"  ' 


but 


and  so   (since   the   inversion  of  the   order   of  integration  and 
differentiation  is  plainly  valid) 
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I 


slnl_  (,^,  ^_  |i2^j^3  =  |(2  ^  ^  cot  1^)  m 


(S     -   T) 


where 
(11.7) 


'(M  =r  <T^^^)'^'  tt!^  ( -'  c  -  ^^>* 


(/; 


'((T 


K 


h/"=  h//'  ^i^/ 


-i<r 


Flgvire   11 


The  integral  (11.7)  is  easily  evaluated,  and  consists  of  two 
parts.   The  first  part  corresponds  to  the  integration  over 
values   s  from  1  to  cos  t,   and  has  opposite  signs  for  K 
on  opposite  sides  of  the  imaginary  axis.   On  the  other  hand 
the  second  part,  which  corresponds  to  the  integration  from 


cos  &     to    oo 


is    single-valued  across    the    iiiiaginary  axis   and 
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so  contributes  nothing  to  the  Integral  (11«6)  because  W" 
and  W'!  are  traversed  in  opposite  directions.   The  contribu- 
tions of  the  Integration  from  1   to   cos  "ii  is 

1   (   3  -  T-  ,1/2   Sin  (^  (^ot/3-lJ:il£)   ds 
J,  ^s  -  cos  t  ^     (g  _  ^j2'     I    3  -  t' 

=  ±  I  3lnp.(   1^  +  I  cot  (3)   r  CSC  (3  .(esc  i^+csc-^^)] 

where  the  plus  sign  is  to  be  taken  on  W"  ,  the  minus  sign  on 

W" .   Hence 
% 

(11.8)  M^=lTra3^sln^-^J^sec  ^.(2  ^  +cot  ^)(^  +  |  cot^). 

[cscp.  (esc  ^  +  CSC  ^  )  ]  dt; 

=  |va^^/   (2csc^^-±^.csc^^L).2^3ec^d!; 


But 


I   (c3c  —^  +   CSC  -^  ).  sec  ^ 


d"t;     =  0 


since  the  integrand  Is  an  odd  function  of  ^  .   Therefore 


(11.9)       M^  =  -anfa'jlij       C3C^5^.  |jsec||d«      . 


"> 
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Finally  by  making  the    substitution  ^(  ^  -  p)   =  z   ,   we   obtain 
(11.10)  My  =  Srr^a^  ^  [     csc^  ^  .    esc  zdz 


w; 


itt 


where     W       is    the  path  of   integration  shovm  in  Pig.  12.      The 
total  virtual  mass      M     is     IL.   -  M      ,   where     M^     la   given  by 


(11.10)    and     li       by   (11.3). 


/,  /t/ 


2L 

2 


Figure   12 
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By   (3«7)    the  dlmenslonless   virtual  mass  Is 


m  = 


Id 


InpA' 


Since  a  =  A  sin  ^  (see  Fig«10),  we  have  from  (11. 3)  and  (11. 10) 
(11.11)  m  =  m^  "  ^z» 


where 
(11.12) 

and 


„  sin'^/a  f    3  (iz 
nu.  =  2  — =r~^   I  C3C  •*=-  .  CSC  zd: 

W 


(11.13) 


2  3 

m  =  — (2  +  3  cos  &  -  cos  /3) 


To  find  C_  (see  (3.9))  It  is  necessary  to  calculate 

P  -a 

the  derivative  of  m  with  respect  to  h  =  x  =  ^  "•"  cos /3  (see 
Pig.  10).  We  have 


(11.14) 


am  _  ^"V       ^"^z 
W  ~  "SB  5b"  * 


where 


(11.15)  ^=-3cot^csc^-m^*i2iin^  .  ^ 


W 


Z   CSC 


//( 


•»^  •  cot^ipcsc  zdz 


and 
(11.16) 


dm 


"55" 


^  =  2  sin^fl  . 
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12.   The  virtual  mass  In  special  cases*   (1)  The  Integral 
(11.12)  can  be  evaluated  in  terms  of  residues  when  the  angle 
P     is  a  rational  multiple  of  tr.   Suppose  that 

(12.1)  yff=l21 

Where     n     and     q     are  relatively  prime   integers   and     n  <  q. 
Then  (11»12)  becomes 

(12.2)  My  ~  ■*  ^^^P  •    ^ 
where 


CSC     —  •  CSC  zd  z 


w'" 


w' 
Let     6      be   a  fijced  positive  nuraber.      Then 

^J     esc'   S£.cscz.  e'^'^dz 

w"'  ^ 

=  -  svim  of  residues   to  the  right  of     W 

=  -     ■  •    sum  of  residues  between  0  and  qtr  • 

Calculating  the   residues  between     0     and     qxr  ,    letting    5   ->  0 

and  observing  that   the  coefficient  of     —    must  vanish,  we  obtain 

6 


(12.3)  I  =  -|"r^(-l)^"^csc3  ^ 


2   ^2 

,(-   1)*'"-'(C3CZ+    ^ 
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If     n  -f  q  Is   odd,   and 

_3 
(12.4)      I  = 


*  ^  2'<-^)'^  [— *  S  6' ^•°^°  ^'L= J 

*='  n  -J 

If  n  +  q  is  even*  Using  these  formalas  a  table  of  values 
for  m  was  computed  by  the  Matheiaatical  Tables  Project 
operating  under  the  Applied  Mathematics  Panel.   These  values 
were  used  in  plotting  the  graph  of  plate  1. 

For  the  purpose  of  checking,  it  is  worth  while  to 
calculate  the  virtual  mass  in  certain  easy  special  cases  where 
it  can  be  calculated  by  other  means. 

(11)  If  ^=|,we  see  from  (12.2)  and  (12.3)  (with 
n  =  1,  q  =  2)  that  m^  =  2  and  from  (11.13)  that  m  =  4  • 
Hence 

This  trivial  case  corresponds  to  the  flow  about  a  sphere  in 
an  Infinite  fluid  (but  only  half  of  the  fluid  is  Involved) 
and  the  virtual  mass  is  then  known  to  be  as  stated. 

(ill)  If  we  let  |3  -»  0  we  obtain  from  ai.l2)  that 


4Tr 


"V"^  *"  -2^ 


°0  4„X 

CSC  zdz  _  .  V~'(-l)'' 


CSC  zdz    J,  \~' 

-73— =  *2. 


^'"   "       ^  =  1    ^ 


3 


f 


W"'  i 


and  from  (11.13) 


m  ^8 
z   3 
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Thus  the  virtual  mass  of  the  fluid  when  the  sphere  is 
tangent  to  the  free  surface  from  below  (  yS  =  0)  is 

00 


-4-^^-1 


r-'       ' 

This  value  agrees  with  the  limit  obtained  when  a  completely 
immersed  sphere  approaches  the  surface  (see  plate  1  and 
Appendix  B),  and  so  is  a  consequence  of  known  formulas  (see  [l]  , 
page  447), 

(iv)  Taking  ySs=  it,  the  virtual  mass  is  obviously  zero. 
However  if  we  keep  a  (instead  of  A)  fixed,  then  we  obtain 
the  mass  for  a  circular  disc  of  radius  a*  Taking  y9  =  ir  in 
(II.  13)  and  (11.10)  we  obtain 

M^  =  2Trpa^.  g^    esc*  zdz  =  ^  p  a^  , 


w'" 


4    3 

and  so  M  =  IL,  -  M  =  "rps.   •   This  is  also  a  known  result 

(see  [1],  page  428). 

13.   Derivative  of  the  virtual  mass  for  rational  ^  •     By  (11.15) 
it  remains  to  evaluate  only  the  integral 

J  =  rr^j-  Z    CSC  ^r   CO*^  ^  CSC  ZdZ  • 

2Wi  j         TT     tr 

w'" 

We  suppose  that    ^=  —^    (n,q)  =  1,   and  (as  above)   let    6    be 
a  fixed  positive  number   so  that 


(13.1)       -L.  f     csc^  ^cot  ^  CSC  ze'^^'dz 
2Tti  J  q  q 


w'" 


i * — -  •sum  of  residues  between  0  and  qtr 

l-(_l)n-»-qe-  <5q^ 
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Calculating  the  residues  between  0  and  qn  and  equating 
the  coefficients  of  o    on  both  sides  of  (13.1),  we  find  that 
If  n  +  q  Is  odd  then 


iLVi-D^'^i   csc^  ilHL  cot  i221 
2/  ^tiw  q       q 

*-  r-' 


+  ^2    Vi'l)^'^  {-LU    esc  2)  +  4  —-.(2  CSCZ)^ 

n 
and  If  n  +  q  is  even 


^  =  ^^2-i^i 


where 


C,  =  TT  r(-i)J-ljcsc=5Ji2i  cot  JHL 


+ 


n 


^2  =  -^  (1  -^  4^  *  1^ XC-D^J^csc^  ifl  cot  i2L 
12n2     qS    2  .^1^  q 


r* 


_^  .2   --    ,,  r  .   o        _2  .3 


q 


,3„.2'<-'^[4(^^csc.,.^^,<A«c.J   ^^ 


z 

n 
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Using  the  fact  that  the  coefficient  of  ^  on  the  right 
side  of  (13.1)  must  vanish,  we  obtain  the  equation 


IT  ' 

(13.2)     T"(-l)^  csc^  J21  cot  iSn. 
jL^  n  q 


-a?    f"(-l)^  FA^CS^  2)  +  4  — ,(C8C  Z)] 


^2 

+ 


n 


V  360  ^4    3g^2    360  / 


This  equation  enables  us  to  evaluate  certain  of  the  sums* 
For  example,  if  in  C,  we  make  the  substitutions  J  =  q  -  j' , 
k  =  n  -  k»,  we  see  that  C,   equals  the  left  side  of  (13.2) 
multiplied  by  the  factor  S-q  • 

Values  for  -r^  used  In  plotting  plate  3  were  computed 
from  these  formulas  by  the  Mathematical  Tables  Project. 
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Appendix  A«      The   Asyrmnetrlcal  Len3» 

14.      Construction  of   the  potential  function*      V/o  construct 
the  velocity  potential  of  the  flow  around  an  asyrametrlcal 
lens  bounded  by  the   spherical   caps    ijj  =a-p  and    ip  =  -y3    , 
where    yS    and    a-/9    are  both  positive.      We  suppose   tliat   the 
lens   Is   immersed  in  an  infinite  fluid  and  that  it  is  moving 
downward   (i.e.    in  the  negative   z-dlrection)   with  unit  velocity. 
The   symmetrical  lens   is   obtained  by  taking   a  =  20    ,   and  this 
is    the   case  which  has  been  discussed  in  the  main  body  of   the 
paper* 


Figiire   13 


The  desired  potential  fvmction    0   is  regular  at  o© 
(  (7=  0,   0  =  0)   and  satisfies    the  boundary  condition 

(14.1)  -|^=-i£ 

on  9n 
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on  (/»  =  a  -  >5  and  (A  =  -  >5  •   The  potential  function  z  satisfies 
(14«1)  but  is  not  regular  at  op  •   Set 

(14.2)  V  =  0-  z  • 

Then     V     satisfies   the   following  conditions: 

(14.3)  -  1^  =  0     on    ^^=a-y?    and       =.p      , 
(14»4)  V     behaves   like   -  z  at  oo  • 

In  order  to  illustrate   a  remark  made   in  section  10, 
namely  that  the   n-fold  Rlemann  space  may  be  used  if   the   inter- 
section angle    (here   equal  to   o«  )   is   rational,   we   suppose   that 
o<  =  ~     where     n     and    m     are  relatively  prime.      The  formulas 
which  we  obtain  are   then  valid  for  all   a.    by  continuity. 

The  potential  function 

(14.5)  V(a.W   =2'['41k«-.)'"'«   -  '*.1„'"'*'] 


-&=■ 


=Z'['2C(k*l)»-.l(''.^'-<''-°i.>.w] 


■fc»l 

has   a  vanishing  normal  derivative   on    i/*  =  -  yS   and    }p=am.p   by 

theorem  (9.1),   and  at  «=(  cr  =  o,  v!*  =  0)   all  terms   in  the   s\wma- 

tlon  are  bounded  except   -q'^'     {o,ip)   -  -  q^'^^      (o,^)   =  -  <^^  {a  xh) 

-2moe  -2niT  o 

which  behaves  at  oo  like    -  z  +  o(l)    (compare    (8.5).   Substituting 
from  (8.2)   into   (l4.5)   and  evaluating  the   sums,   we  obtain  the 
formula 

(14.6)  V(a,W   =^Jq^.   i  [  cot  !!i^j^ -cot  §1  ]    d|  . 


w 
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The    desired  potential  function  is    $  =  V  +  z. 

15.     Virtual  mass   of   the  fluid.      The   calculation  of  the 

virtual  mass    is    similar   to   that  of   section  11.      The  result   is 


(15.1) 

where 
(15.2) 


M  =  My  -  Mj 


=  2na  psin 


U  . 


a 


WT 


CSC    ^  cscz  csc(z  +  ^)dz 


If  3  02  ,         nySv  ^ 

-    ^fJT  °^°     IT  °^°z  C3c(z  -  -~)clz 

W"'    and     W'"    being  the  paths   shown  in  the  diagram,   and 


^/ 


f^. 


f.OC 


Figure   14 
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(15.3)    }A^=  |"Pa^ 


cot^csc^5«(2-co3yS)+  cot-2J^csc^  «9^.(2-co3^a-yS)) 


is  the  mass  of  the  fluid  displaced  by  the  lens. 

If  oe  =  2y3  ,  the  virtual  mass  defined  by  (15.1)  is  of 
course  twice  that  obtained  in  section  11  where  only  half  the 
space  was  considered  to  be  occupied  by  fluid. 
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Appendix  B«   The  Coaipletely  Immersed  Sphere* 


15«   Relation  to  the  lens*   We  consider  the  case  where  the 
sphere  Is  completely  Immersed.   The  sphere  Is  then  moving 
downward  entirely  below  the  siirface  and  does  not  Intersect 
it,  as  shown  in  the  diagram.   We  shall  see  that  the  virtual 
mass  of  the  fluid  in  this  case  is  obtained  from  the  formulas 


of  section  11  by  the  following  sln^jle  rule:  (a) 


M, 


has  the 


<2 


Figure  15 
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same  meaning  as  before;  it  is  the  mass  of  fluid  displaced 

4    "^ 
"by  the  sphere,  therefore  equal  to  -r"^ pk^}    (b)  the  term  My. 

Is  obtained  from  the  formula  (11.9)  merely  by  replacing  the 
real  quantities  f,^and  a  by  1^  ,  1/S  and  la  respectively. 
By  reflection  on  the  plane  z  =  0,  we  see  that  the  flo\» 
corresponds  to  that  around  two  equal  spheres  moving  In  the  line 
of  centers  with  the  same  velocity.   If  the  spheres  are  unequal, 
the  virtual  mass  may  be  obtained  from  the  formulas  of  Appendix 
A  by  the  same  rule.   For  simplicity  we  discuss  only  the  sym- 
metrical case. 

17 .   The  potential  function  and  virtual  mass  of  the  fluid.   S o 
far  a    and  ^  have  always  been  real.   We  now  generalize  by 
allowing  them  to  be  complex.   Write  a  =  A  sin/9,  where  A  is 
a  fixed  positive  number  and  -y9  a  fixed  value  ot   Jp  ,   and  let 
the  complex  bipolar  coordinates  a ,  ip    "be   connected  with  the  com- 
plex coordinates  r,z  by  the  relations  (compare  (6.3)) 

^       ^  sinh  cr  ^       „  sin  ^ 

r  =  a     "  ,      z  =  a         '  • 

cosh<T  -  cos  ip  cosh  a  -  cos  ip 

Let    6    be  real  and  set     x  =  r   cos  0     ,   y  =  r  sine  ;    then  x,y,z 
are   complex   cartesian  coordinates   in  space  which  becomes   real 
when    o  ,^/^    and   yS    are  all  real  or  all  pure   imaginary.      In   the 
flrs^  case   o  ^  xp ,   6   are   toroidal  coordinates,    in  the   second 
Im  (7    f   Im  ip   ,    6     are  the  so-called  bipolar  coordinates   in  space. 
The   locus     xp  =  - 0      is   a   spherical    cap  in  the   first   system  and 
a   sphere   in  the   second;   in  both  systems    the  radius   is     A     and 
the   center   is   at     r  =  0,    z  =  -  A  cos  yS    • 

The  function     q^  (o,ip)   defined  by   (7.3)   remains   a  poten- 
tial function  even  though    o    and     ^  are   complex.      However  we 
are   interested  in  this   generalization  only  in  the   case  where 
a   ,  ip  f   ^Q     ^^^  fi    ^r®  ^11  pvire   imaginary,    a   condition  we  suppose 
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Is  satisfied  henceforth.   In  order  to  have  r  S  0  we  suppose 
also  that  cr  =  l|cr(  where  lo-l  <  v    • 

Since  a    and  Y  are  imaginary  all  the  branch  points  of 
q^(«r,'S^  (as  a  function  of  ^  )  lie  on  the  line  Im  ^  =  Im-v^.  For 
each  Integer  k,  k  =  0,  ±1,  ±2,  •••,   we  cut  the  fe  -plane 
along  the  straight-line  segment  1.   connecting  the  two  adja- 
cent branch  points  y-  icr+  Skrr  and  y  +  1  (r+  2(k+l)Tr  ;    in 
the  plane  so  cut  q^   is  single-valued.   Suppose  first  that 
|cr|  >  0  and  let  L.  denote  a  closed  contour  surrounding  1. 
as  in  the  diagram.   We  write  W  =^^1'.  . 


^-a                           Z_, 

r    ^.          ^. 

.—  )        (, .)    . 

"C- 3     c  — 

yr-ir-sn           t/f^itr 

Vr-iT               Vr^ia--t-ZV 

Figure  16 

We  suppose  that  the  sphere  is  moving  downward  (negative 
z-direction)  with  unit  velocity.   Then  the  desired  potential 
function  $  is  regular  at  <o   and  satisfies  the  conditions  that 
^  =  0  on  '>|^=  0  and  "  ^  =  "  "35  on  *^®  sphere  '^=  -/3  • 
If  we  set  V  =  $  -  z,  it  is  easy  to  verify  that 


(17.1) 


V  = 


-otJ 


ff  °^^^ 


w 
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Since  p  is  imaginary,  we  observe  that  the  integrand  tends 
exponentially  to  zero  as   fe,->  °o   on  W. 

The  calculation  of  the  virtual  mass  M  of  the  fluid 
is  formally  similar  to  that  of  section  11.   The  final  result 
is  that 

(17.2)  M  =  M^  -  M^ 
where 

(17.3)  ^z'^J^^^^ 
and 

(17.4)  H^  =  .  2Trea'  ^  J  csc^  ^'^   sec  ?|  at  , 

the  integration  being  along  the  real  axis  in  the  negative 
direction.   Since  the  integration  in  (11.9)  is  along  the 
imaginary  axis  in  the  positive  direction,  we  obtain  the  rule 
stated  in  section  16. 

Finally,  calculating  the  residues  in  the  upper  half -plane 
we  obtain  for  the  dimensionless  virtual  mass  m  =  M  /(iw^A^) 
the  formula 

oo  o    -I 

(17.5)  m  =  I  *  4  V-(,i)kr   sin  (3 

-f^      Lsln(k  +  a)p 

(Here  /3  is  of  course  imaginary. ) 

As  (3  varies  from  ir  to  0  along  the  real  axis  and 
then  up  from  0  to  ioo  along  the  imaginary  axis,  the  sphere 
passes  through  the  surface  and  then  moves  down  through  the 
fluid  to  infinite  depth.   The  virtual  mass  m  and  its  deriva- 
tive ~  vary  continuously,  but  we  see  from  plate  2  that 
^-Bi  has  a  discontinuity  at  |3  =  o  (b  =  1  +  cos  p  =  2). 
db^ 
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